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OPTIMAL  DESIGN  OF  LOCALLY  ORTHOTROPIC 
ELASTIC  FLAT  BODIES  WITH  WEAK  BINDING 

G.  I.  Bryzgalin 
(Volgograd) 

Solved  are  problems  of  optimal  design  for  elastic  bodies  made 
from  a  material  such  as  fiberglass  reinforced  plastics  or  metal 
reinforced  by  fibers  with  a  high  elasticity  modulus  [1,  2],  It  is 
assumed  that  the  composite  components  are  in  a  plane  strain  or  a 
plane  stress  state  and  the  reinforcement  is  directed  along  two 
mutually  orthogonal  families  of  plane  curves.  The  equations  for 
these  curves,  the  reinforcement  strength  at  each  point,  and  the 
body  thickness  are  determined  during  the  designing  for  a  given 
boundary  profile  and  conditions  on  it. 

Let  us  present  a  rather  small  body  element  cut  by  planes 
normal  to  the  directions  of  reinforcement.  On  the  two  opposite 
edges  orthogonal  to  the  first  direction  of  reinforcement  we  will 
apply  expanding  or  compressing  normal  stresses  a^.  We  will  assume 
that  the  material  is  deformed  as  a  monolith,  i.e.,  deformations 
along  the  stress,  on  the  average,  are  equal  at  all  points.  Let 

be  the  stress  in  the  reinforcement  fibers  of  the  first  direction 
s-^  -  cross  sectional  area  of  these  fibers  per  unit  area  of  material 
and  E*  -  elasticity  modulus  of  fibers. 
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If  the  edge  of  the  examined  element  has  a  unit  area,  then  the 
remaining  part  of  material,  composed  of  binder  layers  and  fibers 
oriented  in  a  transverse  direction,  will  occupy  on  it  area  1  -  s^. 
The  average  elasticity  modulus  of  this  remaining  part  is  designated 
in  terms  of  E",  and  the  average  stress  -  Using  the  well-known 

formula  of  mechanical  mixing  we  obtain 


t  o»"(l  «—#»).  i>/  «  E't\,  or  i-*  E"rt 


CD 


Considering  that  s.^  does  not  equal  zero  and  E"  is  considerably 
less  than  E1,  we  can  write  the  approximate  equality 


H|  »  o/l|  *  E'itft 


(2) 


Calculations  based  on  more  specific  formulas  [3]  indicate  that 
if  the  modulus  of  elasticity  of  binder  Ec  is  ten  times  smaller  than 
the  reinforcement  modulus,  then  the  replacement  of  equality  (1)  with 
equality  (2),  at  worst,  yields  an  error  (contributing  to  the  safety 
factor)  on  the  order  of  30%  of  the  maximally  possible  stress  in 
the  material. 


The  following  experimental  factors  also  speak  in  favor  of  this 
substitution: 


a)  the  elasticity  modulus  of  the  binder  of  fiberglass  rein¬ 
forced  plastics  is  10-20  times  lower  than  the  elasticity  modulus 
of  reinforcement; 

b)  usually,  binders  are  subject  to  considerable  creep,  and 
coefficient  E"  in  the  preceding  formulas,  in  essence,  is  not  a 
constant  value  but  an  operator,  such  that  with  =  const  11m  E"e^  = 
=  0  when  t  -*■<*>  (t  —  time).  So  that  for  parts  working  under  condi¬ 
tions  of  prolonged  action  of  slightly  changing  stresses,  relation¬ 
ship  (2)  can  prove  to  be  more  acceptable  than  (1); 

c)  the  range  of  the  elastic  and  strength  curves  of  composite 
materials  is  tens  of  percent  from  sample  to  sample. 
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The  elasticity  modulus  of  such  an  element  can  be  considered 
as  identical  during  expansion  and  compression  [4],  and  Poisson 
coefficients  during  deformation  along  the  fibers  as  small  (so  small 
that  they  can  be  disregarded);  therefore,  we  extend  the  relation¬ 
ship  of  type  (2)  to  the  general  case  of  plane  stressed  state. 

We  will  limit  ourselves  to  the  designs  in  which  the  main  axes 
of  stress  and  strain  tensors  coincide  with  the  fiber  directions  of 
reinforcement  at  every  composite  point  (rational  designs  [2]). 

Then  in  the  curvilinear  orthogonal  system  of  coordinates  a3  along 
whose  coordinate  lines  the  reinforcing  fibers  are  placed,  the 
following  law  governing  the  connection  between  strain  and  stresses 
can  be  written  as: 

Oi  n  a-  *=  £«. f.  ( 3  ) 

Here  o^,  a 2  -  principle  stresses,  e^,  e2  —  principle  strains, 

E  —  constant  selected  on  the  basis  of  experiments  on  the  composite 
material. 

The  problem  concerning  the  strength  criterion  of  such  hypo¬ 
thetical  material,  in  other  words,  concerning  the  limits  of  appli¬ 
cability  of  relationships  (3)  is  rather  complex  and  it  is  doubtful 
whether  it  can  have  a  unique  answer  suitable  for  description  of 
all  actual  materials  of  a  similar  type. 

Taking  advantage  of  the  fact  that  tensile  strengths  of  com¬ 
posite  materials  during  extension  and  compression  are  values  of 
the  same  order,  we  will  assume  that  the  permissible  state  of  an 
element  is  characterized  by  deformation  values  not  exceeding  a 
certain  maximum  value  of 

M«S*.  I»:i  ^  t  (^) 

For  reinforcement  which  remains  elastic  right  up  to  destruction 
(oj  =  E,ei,  i  =  1,  2),  this  means  that  the  maximum  value  of  internal 
stress  in  fibers  determines  the  material  strength. 
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Of  course,  during  compression  or  during  a  complex  stressed 
state,  destruction  occurs  not  necessarily  due  to  a  break  in  fibers; 
however,  criterion  (4)  can  be  used  safely  if  the  materials  tested 
have  different  s^,  s2,  if  we  construct  families  of  strength  curves 
at  different  stressed  states  in  plane  e^e2,  and  se^ect  the  square 
j.  of  the  permissible  states  such  that  none  of  the 

strength  curves  should  intersect  it  (Fig.  1). 


~st' 


s'. 


Fig.  1. 


Since  there  are  no  fibers  in  a  direction  normal 
to  the  coordinate  lines  of  system  ct3,  then  within 
the  accepted  permissible  limits  we  can  assume  that 
the  third  principal  stress  is  equal  to  zero  c^  =  0. 


The  stress  components  with  mixed  indexes  equal  zero  due  to  the  fact 
that  this  is  precisely  how  the  rational  designs  are  constructed; 
so  that  directions  a  -  const  and  3  =  const  were  the  principal 
tensions  of  stress  and  strain  tensors  o12  =  =  ao 3  = 


In  view  of  this  fact  there  will  not  be  any  fundamental  differ¬ 
ence  between  the  plane  stress  and  plane  strain  states.  We  will 
only  note  that  with  a  plane  strain  state  the  design  is  for  a  com¬ 
posite  layer  with  a  unitary  thickness,  while  with  a  plane  stress 
state  -  the  thickness  can  be  arbitrary  and  even  vary  from  point  to 
point;  thus,  only  its  minimally  possible  value  is  selected  during 
designing.  The  problems  will  become  wholly  identical,  if  in  the 
latter  case  in  place  of  stresses  and  levels  of  reinforcement  we 
introduce  values 


(h  —  composite  thickness  at  a  given  point) 


(5) 


If  we  consider  that  the  above  examined  volume  element  has 
}  ’  ght  h  and  unitary  transverse  dimensions,  then  T^,  T2  -  forces 
acting  on  this  element  along  directions  a  and  3;  S^,  S2  —  fiber 
volumes  of  corresponding  direction  in  it.  (For  simplicity,  T^, 


T2  will  also  be  called  stresses.) 


The  equations  of  equilibrium  of  such  volume  element  ([5], 
Chapter  2)  can  be  obtained  from  (8.2),  if  we  assume  that,  there. 


°12  '  0 


*  _  oB-.  d  oUs 

— 'iH'.Ti)—  Tz  •  0.  — (HJ «)—  —T%  ■  0 

9u  UU  /*J»  nf} 

(H^,  H2  -  Lam£  parameters  of  system  a0) 


(6) 


Formula  (8.2)  in  work  [6]  (Chapter  2),  when  e^2  =  0,  leads  to 
the  following  conditions  of  deformation  consistency: 


Equations  (6),  (7)  should  include  boundary  conditions  on  body 
profile,  conditions  (4),  and  positive  conditions  of  reinforcement 
levels  >  0,  S2  >.  0. 

Let  the  problem  be  solved  with  regard  to  the  theory  of  elas-  ’ 
ticity  for  an  isotropic  body  whose  middle  plane  occupies  area  n 
bound  by  contour  L  on  plane  a0.  Then  for  constructing  a  reinforced 
body  design  with  the  same  boundary  profile  and  same  conditions  on 
it,  we  can  assume  that  the  reinforcement  fibers  are  directed  along 
the  lines  of  principal  stresses  of  the  isotropic  body  (theorem  on 
the  existence  of  rational  design  [2]).  For  the  selected  system  of 
coordinates  aB  (whose  coordinate  lines  coincide  with  the  isostatic 
grid  of  the  isotropic  body)  we  have  nine  equations  (3),  (5)-(7) 
and  ten  functions 

Tt.  f;.  O I.  <1;  Si.  S;.  «|.  t:  >1.  I;  (  8  ) 


The  initial  form  of  these  functions  and  the  body  thickness  are 
determined  according  to  the  indicated  theorem,  and  beyond  this  it 
is  possible  to  state  a  variation  problem,  i.e.,  require  that 
functions  (8)  yield  an  extremum  to  a  certain  functional.  In  this 
work,  the  design  withstanding  given  loads  and  having  the  smallest 


volume  of  reinforcement  among  all  rational  designs  we  will  consider 
as  optimum.  In  other  words,  for  the  optimum  design,  functions  (8) 
satisfy  all  the  conditions  written  earlier  and  yield  a  minimum  for 
the  following  functional: 

f<S,  +  S:i<IQ  -  mit  (Q) 


In  a  specific  case,  when  both  principal  stresses  in  an  iso¬ 
tropic  body  have  the  same  sign,  the  design  is  constructed  on  the 
basis  of  the  corresponding  theorem  of  existence,  in  this  case 

•  l  “  *1  “  1*  “  fu#«t  (1C 


Another  specific  case  corresponding  to  the  following  conditions: 


/»  »  -l:  »  »  ron«t 


is  indicated  in  work  [7],  The  possibility  of  constructing  a  design 
which  would  satisfy  equalities  (11)  is  restricted  by  the  condition 
of  deformation  consistency.  Actually,  assuming  that  =  const, 
e2  *  const  in  (6)  and  using  the  condition  of  equality  to  zero  of  the 
tensor  of  curvature  of  space  <x8  in  the.  initial  state,  we  find 


H  ii 

</<(z  lUg  'III 


i»&*  •••;  ft* 


In  this  case  the  second  equality  proves  to  be  a  symmetrical 
substitution  of  indexes  and  coordinates;  furthermore,  from  it  we 
can  easily  obtain 

w/./cf. -/,«)//,  (13) 

where  f(a)  —  arbitrary  function.  Using  the  known  expressions  for 
the  curvatures  of  families  of  plane  curves 


H,H,  <><i 


we  rewrite  (13)  thus 


— x.ff.  —  /(o) 


Condition  (14)  was  shown  in  [7];  geometrically  it  mpans  that 
the  increase  in  the  inclination' angle  of  a  tangent  to  line  a  per 
single  increase  of  coordinate  a  does  not  depend  on  0.  '  i 

Generally,  we  have  to  solve  the  variation  problem;  moreover, 
we  should  assume  that  functions  T2,  S2  can  have  a  break  along  line1 
a  =  const,  and  functions  ,  S ^  -  along  line  8  =  const.  After  we 
determine  functions  $1,  S2,  thickness  h  can  be' found  in  the 
following  manner.  Let  !  1 

.  niux  (5*(u.  1)  -r  $:C«.  p)  H  V  »  , 

O 

\ 

At  point  (aQ,  8q),  where  this  maximum  is  achieved1,  we  will 
assume  that  s^  +  s2  =  s*  (s#  —  greatest  possible,  reinforcement  1 
density,  s-  <.  1).  Then 

*  v 

^  *  A  / u,  **  **  •  Si/  .V.  $t  ■  /  A  (15) 

1  . 

There  are  other  possibilities  for  selecting  the  thickness  and 
level  of  reinforcement'  since  addition  or  Removal  of  pxcess  binder 
does  not  (within  the  framework  of,  this  system),  affect  the  carrying 
ability  of  the  design.  _  , 

i  ,  , 

i  r 

Further,  we  are  presenting  examples  illustrating  various  possi 

bilities  of  spatial  optimum  designs.  1 

> 

A.  Half-plane  loaded  on  the  edge  section  by  a'  normal  uniform 
force  of  intensity  P.  First  of  all  let  us  investigate  a  stressed 
state  in  the  corresponding  is'otropic  body.  From' formulas  (4,  4’, 

5)  of  work  [8],  (page  352)  we  can  obtain  the  following  expression 
for  the  tensor  components  of  stresses  in  an  elastic , isotropic 
half-plane  loaded  in  the  same  way  (there  is  an  error  in  the  corre¬ 
sponding  equalities  [8]) 

«  .t'/I  Q..|/r,rii  . 

T,  —  o.*  n  —ill,— dil— or>  (0,  +  0:.  /  f,f»J 
t,.  -  n„li  ™  (8,  +  »l;|  r,r, 

(notations  are  in  Fig.  2).  in  this  case  it  turns  out  that  the 


(16) 


I 


Figure  2. 


principal  direction  a  forms  angle  y  =  1/2 ( © ^  +  axis  x,  and 

the  principal  stresses  are 

T,  .t'PI— +  'r.r.J.  r,  -•  n-*P|— ft  —  ‘lay /rtt.)  (  17  ) 

Here 

ft  —  81  —  61.  8,  —  «rctf  1— yl  <»-«)!.  U,  -  tn  ig  (—»/<*  +  a)]  ( 1 8  ) 

;  Expression  (17)  can  be  presented  in  the  form 

T,  —  n-'P[— 8  —  tiaft],  f«  —  l— ft  +  »in  ftl  (19) 


.  It  is  obvious  from  (19)  that  both  principal  stresses  are 
always  negative.  Consequently,  it  is  possible  to  construct  a 
design  with  a  reinforcement  which  is  uniformly  stressed  and  with 
identical  principal  strains 

, 

i«  »,  •  (»’»*!  (20) 

The  trajectories  of  principal  stresses,  along  which  the  rein¬ 
forcement  fibers  are  oriented,  in  this  case,  are  families  of  co- 
focal  ellipses  and  hyperbolae  with  foci  at  points 

/  «  n  *  *  —  a.  ».  <* 

here,  the  hyperbolae  determine  direction  a. 


On  the  x-axis  segment  [-a,  a]  the  principal  stresses  are 
maximum,  under  condition  (20)  this  indicates  the  maximality  of 
value  S1  +  P-, ,  then,  according  to  (15),  we  find  that 
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It  is  interesting  to  note  that  if  the  load  on  segment  [-a.  a] 
is  not  constant  and  the  constant  displacement  v^  =  const  (smooth 
rigid  stamp)  with  the  same  total  force  2Pa,  then  the  optimum  design 
coincides  with  the  constructed.  Actually,  the  found  deformation 
field  permits  us  to  satisfy  condition 

i  (>,  i'|  **  «•'  (»,  —  »!,<•>*«  —  t.  -j-  r'| 

It  is  evident  from  (19),  (21)  that  when  r^,  r0  -*•  «,  s^,  s2 
0,  consequently,  at  the  sufficient  distance  from  the  loaded 
segment  the  reinforcement  can  be  terminated  (when  the  contribution 
of  depleted  reinforcement  to  the  total  material  strength  proves  to 
be  small  in  comparison  with  the  actual  strength  of  the  binder). 

B.  Circular  plate  loaded  with  a  uniform  tangential  force  along 
the  internal  and  external  contour .  We  construct  a  design  which 
satisfies  conditions  (11).  As  a  matter  of  course  we  introduce 
polar  coordinates  rO  and  accept  for  the  reinforcement  directions 
condition 

V  *=-  0  x  */»n 

(plus  for  a  and  minus  for  3) 

specifying  that  direction  a  forms  angle  1/4tt  with  the  radius-vector, 
and  angle  y  with  axis  x.  In  this  case  the  equation  of  line  a  in 
system  xy  has  the  form 

y"  nr  J,  —  ■>)  /  (j  -  •  *  | 

A  solution  of  the  latter,  written  in  polar  coordinates  6  =  In  cr 
(logarithmic  spiral),  is  c  =  const.  Equation  of  lines  3:  6  = 

=  -In  cr. 

Families  of  such  curves  s££A»£y  •ewidTt'ian  (14).  The  tensor 
components  of  stresses 

r„  -  /v  /  r».  r.  -  r, « 

Here  F  -  force  per  unit  length  of  the  internal  contour,  a  - 
internal  plate  radius,  the  external  radius  can  be  arbitrary. 
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The  principal  stresses,  thickness,  and  reinforcement  levels 
are  expressed  by  formulas 

From  conditions  (11)  we  can  also  determine  the  displacements 
along  the  radius  and  in  the  circumferential  direction 

■'  —  0,  u*  -  2»r In  (f  /•) 


C.  Circular  plate  under  internal  pressure .  The  radial  and 
circumferential  directions  we  take  to  be  the  reinforcement  direc¬ 
tions,  then  the  equations  of  equilibrium  and  compatibility  will 
assume  the  form 

r,  —  <rP)*.  »,  -Iff,)’  (22) 

Here  subscript  1  pertains  to  the  radial  conditions  and  sub¬ 
script  2  —  to  the  circumferential;  the  prime  indicates  differentia¬ 
tion  with  respect  to  r. 


It  is  shown  in  [1]  that  in  the  absence  of  external  pressure 
the  construction  of  a  design  with  e-^  ®  e2  *  const  is  impossible 
according  to  (10).  It  is  readily  seen  that  construction  of  a 
design  which  satisfies  conditions  (11)  is  also  impossible  due  to 
the  fact  that  the  second  condition  (22)  is  not  satisfied. 


We  will  assume  that  the  reinforcement  fibers  are  oriented, 
this  leads  to  the  following  condition  concerning  the  reinforcement 
level  in  radial  direction  [1]: 


£in|i/r  (m  —  eon»H 


Boundary  conditions 


r,  ~  -P  when  r  -*  t.  T,  «■  »  wicn  •  k 


Expressing  all  variables  Ifi  "berms  o?  e2,  we  find  specifically 


.V  .'1  /  «: 
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We  select  S2  such  that  volume  (9)  of  reinforcement  is  the 


smallest,  i.e.. 


.  r»n"  +  -rti 

\ -  — .  dr  e=  mm 

J  fj 


The  Euler  variation  equation  has  the  form 


rt,t,"  +  —  0 


Its  general  solution  is  £2  =  ci  exP  where  =  0,  3^^  ^ 

^  0  when  r  =  b;  consequently  (see  (3),  (5))>  we  have  e-^b)  =  0. 

From  this  c2  =  b.  With  aid  of  (22)  we  obtain 

»»  -»  f,np(6/r)  ti  ■»  t,U  —  b/r)  exp  (t/r)  (26) 

Let  b  <.  2a  (for  b  >  2a  the  investigation  is  similar  but  the 
expression  for  constant  c^  is  different). 

Functions  (26)  have  the  greatest  values  when  r  =  a,  here 

mkS<M.  If;!!  - 

We  will  assume  that  this  maximum  is  equal  to  the  maximum 
deformation  value  e  and  thus  determine  the  constant 

r<  =  rnp(—t/t) 

Now,  from  (2*1)  we  determine  the  level  of  reinforcement  S2(r), 
which  yields  the  extremum  to  the  examined  functional 

S:(r)  -  Pk;/r* 

The  remaining  boundary  condition  we  use  to  find  constant  p 
finally. 

Si  -  Pa * Et  <»  —  «)r,  S,  <-  Pa’-f- 
*-*/*(•*  +  »»>/£«.(»-«)• 

The  design  is  constructed.  We  can  show  that  the  second 
variation  of  functional  (25)  is  an  extremely  positive  functional 
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at  "point”  c2  ■  o i  exp  (b/r);  consequently,  the  found  extremum  io 
a  minimum. 

For  reinforcement  volume  (9)  we  obtain  the  following  formula: 

V*  -  2a*4(«  + »)  /«• 

In  conclusion,  we  will  examine  the  behavior  of  representative 
point  (e^  c2)  describing  the  stressed  state  of  the  body  elements 
along  a  certain  plate  radius.  For  this  we  will  express  in  terms 
of  e2(b  -  2a)  with  the  aid  of  (26) 

•l  -  -l|(l +  !■<*! /«)) 


A  corresponding  curve  is  shown  in  Fig.  2.  When  r  changes  from 
a  to  2a  the  representing  point  travels  the  curve  from  B  to  Q. 

As  we  can  see,  there  is  no  analogy  here  with  the  optimum 
design  of  the  rigidly  elastic  bodies,  where  it  is  most  convenient 
and  always  possible  to  use  the  points  on  a  restricting  surface. 

In  the  examined  case,  the  points  at  which  the  limiting  curve  inter¬ 
sects  the  axes  of  coordinates  (for  example)  are  less  suitable  than 
many  internal  points,  since  one  of  the  reinforcement  directions 
for  them  is  not  loaded. 

The  author  expresses  his  appreciation  to  Yu.  N,  Rabotnova 
for  his  valuable  suggestions. 
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